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The Tevatron experiments DO and CDF have found evidence for single top quark production, 
based on datasets between 0.9 fb _1 and 2.2 fb _1 . Several different multivariate techniques are 
used to extract the single top quark signal out of the large backgrounds. The cross section 
measurements are also used to provide the first direct measurement of the CKM matrix element 
|V«|- 



1 Introduction 

Evidence for single top quark production at the Tevatron and a first direct measurement of 
the CKM matrix element \Vtb\ was first reported by the DO collaboration 0} In cont rast to top 
quark pair production through the strong interaction, which was observed in 1995^^1, single 
top quarks are produced via the weak interaction. The Feynman diagrams for standard model 
(SM) s-channel (tb) and t-channel (tqb) single top quark production are shown in Fig. [TJ There 
is third production mode, associated production of a top quark and a W boson, but its cross 
section is so small that it will not be considered further. The SM cross s ectio n for the s-channel 
process pp^tb + X, tb + X is 0.88 ± 0.14 pb at NLO for m top = 175 GeV^. At the same order 



and mass, the cross section for the t-channel process pp—>tqb + X, tqb + X is 1.98 ± 0.30 pbSEl. 

Measuring the single top quark production cross section provides a direct measurement of 
the CKM matrix element | Vtb \ • The single top quark final state also allows for studies of the top 
quark polarization, and it is sensitive to many models of new physics, for example flavor changing 
neutral currents via the gluon^ or heavy new bosons W that only couple to quarks^. The s- 
channel process is also an important background to Higgs searches in the associated production 
mode, and the advanced analysis techniques used in the single top searches will be applicable 
to Higgs searches as well. 
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Figure 1: Feynman diagrams for s-channel (left) and t-channel (right) single top quark production at the Tevatron. 



The DO collaboration has updated two of their analysis methods using a dataset of 0.9 fb . 
The updated results, including a combination of different methods are presented below. The 
CDF collaboration has analyzed a dataset of 2.2 fb -1 and significantly improved the sensitivity 
to single top quark production. These new results are presented below. 

2 DO results 

2.1 Event selection 

The DO analysis selects electron+jets and muon+jets events in 0.9 fb -1 of data with the following 
requirements: One high-p^ lepton (electron (pT > 15 GeV) or muon (pr > 18 GeV)), missing 
transverse energy Et > 15GeV, and between two and four jets with jet pr > 15 GeV and 
jet 1 pt > 25GeV, at least one is tagged with a neural-network based b-tagging algorithm. 
Additional cuts remove fake-lepton background events. Events are collected by lepton+jets 
trigger requirements. 

The number of events observed in data and expected from the background model and SM 
signal is shown in Table [TJ The largest sources of systematic uncertainty are the background 
normalization, jet energy scale, as well as b-tag and trigger modelling. 

Table 1: Numbers of events expected by DO in 0.9 fb - for electron and muon, 1 fo-tag and 2 fe-tag channels 

combined. 





2 jets 


3 jets 


4 jets 


s-channel 


16±3 


7±2 


2±1 


t-channel 


20±4 


12±3 


4±1 


tt 


59±14 


134±32 


155±36 


W+jets 


531±129 


248±64 


70±20 


Multijets 


96±19 


77±15 


29±6 


Total background 


686±131 


460±75 


253±42 


Data 


697 


455 


246 



Table [T] shows that after selection cuts, the expected SM single top signal is small compared 
to the background sum, and in fact the signal is significantly smaller than the background 
uncertainty. Thus, more advanced techniques are required to extract the signal. 

2.2 Multivariate techniques 

The DO analysis employs three different multivariate techniques to extract the single top quark 

signal out of the large backgrounds. The boosted decision tree (BDT) analysis has not changed 

[7l 

since the publication of evidence for single top quark production . Here we focus on the Bayesian 
neutral network analysis and the matrix element analysis, both of which have been re-optimized. 
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Figure 2: Comparison between data and background sum for the Bayesian neural network output. Shown is 
the full distribution (left), and the high-discriminant region (right). The signal has been normalized to the SM 

expectation. 



In a conventional neural network, the network parameters and weights are determined in an 
optimization (training) procedure. Rather than optimizing for these network parameters once 
and then fixing them, the optimal network configuration can be obtained as an average over 
many different values for the network parameters. In this Bayesian procedure, an integration 
over all of the possible network parameter space is performed. The network architecture is fixed, 
and the weight of each set of parameters is obtained through a Bayesian integration. The final 
network discriminant is then the weighted average over all the individual networks. Fig. [2] shows 
the output of the BNN for the DO data. 

The Matrix element analysis starts from the Feynman diagrams for the single top quark 
processes and uses transfer functions to relate the parton level quark-level information to the 
reconstructed jet and other information. Matrix elements for the single top quark signal as 
well as the VF+jets backgrounds are included. For 3-jet events, a top pair matrix element is 
also included. For each event, an integration over the phase space is performed, employing the 
transfer functions to compute the probability for this particular event to arise from a specific 
matrix element. A likelihood function is then formed as the ratio of the signal and signal plus 
background probabilities. 

2.3 DO summary 

The cross section is measured as the peak of the Bayesian posterior probability density, shown 
in Fig. [3] for the ME analysis. The three different methods measure the following cross sections 
for the sum of s- and t-channel: 

a obs ( pp _^ tb + x ^ tqb + = 4.9+1-4 pb 

= 4At\i pb (BNN) 
= 4.81} j pb (ME). 

The measured cross sections are consistent with each other and above the SM expectation. 
The decision tree analysis has also measured the s- and t-channel cross sections separately, 

o- obs (pp -> tb + X) = 1.0±0.9pb 
o-obs _^ tqb + = P b > 

where the standard model cross section is used for the single top process not being measured. 

Removing the constraint of the standard model ratio allows to form the posterior probability 
density as a function of both the tb and tqb cross sections. This model-independent posterior is 
shown in Fig. [3] (right)for the DT analysis, using the tb+tqb discriminant. The most probable 
value corresponds to cross sections of a(tb) = 0.9 pb and a{tqb) = 3.8 pb. Also shown are 
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Figure 3: Posterior probability density for the matrix element analysis as a function of the sum of s-channel and 
t-channel cross sections (left), and for the BDT analysis as a function of both the s-channel and t-channel cross 

sections (right). 



the one, two, and three standard deviation contours. While this result favors a higher value for 
the f-channel contribution than the SM expectation, the difference is not statistically significant. 
Several models of new physics that are also consistent with this result are shown in RefM These 
updated results have recently been published™. 



3 CDF results 

3.1 Event selection 

The CDF analysis selects electron+jets and muon+jets events in 2.2 fb _1 of data with the 
following requirements: One high-pr lepton (px > 20 GeV), E-t > 25GeV, and two or three jets 
with jet px > 20 GeV, at least one of which is tagged by a displaced vertex tagging algorithm. 
Additional cuts remove fake-lepton background events. Events are collected by single- lepton 
trigger requirements. The matrix element analysis uses additional triggers in the muon channel 
to increase the acceptance. 

The number of events observed in data and expected from the background model and SM 
signal is shown in Table [2j The largest sources of systematic uncertainty are the background 
normalization, jet energy scale, and b-tag modelling. Again, it is clear that a advanced analysis 



Table 2: Numbers of events expected by CDF in 2.2 fb 1 for electron and muon, 1 6-tag and 2 6-tag channels 

combined. 





2 jets 


3 jets 


s-channel 


41±6 


14±2 


t-channel 


62±9 


18±3 


tt 


146±21 


339±48 


V7+bottom 


462±139 


141±43 


W^+charm 


395±122 


109±34 


VF+light 


340±56 


102±17 


Z+jets 


27±4 


11±2 


diboson 


63±6 


22±2 


Multijets 


60±24 


21±9 


Total background 


1492±269 


755±91 


Data 


1535 


752 



techniques are required to extract the signal. 
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Figure 4: Comparison between data and background sum for the t-channel likelihood discriminant (left), the 
neural network discriminant (center), and the light quark jet pseudorapidity in the high-discriminant region for 
the neural network analysis (right). The signal has been normalized to the SM expectation. 

3.2 CDF Likelihood Function 

A multivariate likelihood is built from several kinematic variables that each separate the single 
top quark signal from the backgrounds. One special variable is a specially developed b-tagging 
neural network that aids in separating b-quark jets from light quark and c-quark jets. An 
additional special variable is a kinematic solver using constraints from the W boson mass and 
the top quark mass to determine if an event is well reconstructed. Another special variable is 
the t-channel matrix element, which uses the kinematic information provided by the kinematic 
solver. The likelihood discriminant for the t-channel likelihood is shown in Fig. H] (left). 

The measured cross section is obtained as the peak of a Bayesian posterior probability. The 
likelihood analysis measures a cross section of a(tb+tqb) = 1.8igJ pb, below the SM expectation. 

3.3 CDF Neural Network 

Several kinematic variables as well as the b-tagging neural network output are combined in a 
neural network. Four different networks are built with 10-14 variables each, trained separately 
for 2-jet and 3-jet as well as 1-tag and 2-tag events. The full neural network output distribution 
is shown in Fig. [J] (center), and the signal region is shown in Fig. [J] (right). The neural network 
analysis measures a cross section of a(tb + tqb) = 2.0]t 'g pb, below the SM expectation but 
consistent with the SM within uncertainties. 



3.4 CDF Matrix Element 

The matrix element analysis uses the same approach as described above, but also includes a top 
pair matrix element in the 2-jet bin. The matrix element for top quark pair events has more final 
state particles than the single top process, and these additional particles have to be integrated 
out. This is done by integrating over the kinematics of the hadronically decaying VF-boson in a 
lepton+jets top pair event. 

The Bayesian posterior probability density for the Matrix element analysis is shown in Fig.El 
showing the measured cross section and the measurement uncertainty. The mesured cross section 
is o~(tb + tqb) = 2.2±°;^ pb, again below the SM expectation but consistent with the SM within 
uncertainties. The CKM matrix element \Vtb\ is also extracted from the posterior probability 
and a lower limit is found to be \V t b\ > 0.59 at the 95% confidence level. 
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Figure 5: Data-background comparison for the matrix element discriminant (left) and Bayesian posterior density 

distribution observed by the Matrix element analysis. 



4 Summary 

Both Tevatron experiments have found better than 3 sigma evidence for single top quark pro- 
duction and have made the first direct measurement of the CKM matrix element \Vtb\ using 
advanced multivariate techniques. The CKM matrix element |V#>| can be measured to better 
than 15%. Further improvements to the analyses are in progress and both experiments are 
working towards observation of single top quark production at the 5 sigma level. 

Acknowledgments 

We thank the Fermilab staff and the technical staffs of the participating institutions for their 
vital contributions. 



References 

1. V.M. Abazov et al. (DO Collaboration), Phys. Rev. Lett. 98, 181802 (2007). 

2. F. Abe et al. (CDF Collaboration), Phys. Rev. Lett. 74, 2626 (1995). 

3. S. Abachi et al. (DO Collaboration), Phys. Rev. Lett. 74, 2632 (1995). 

4. Z. Sullivan, Phys. Rev. D 70, 114012 (2004). 

5. Q. H. Cao, R. Schwienhorst and C. P. Yuan, Phys. Rev. D 71, 054023 (2005) 
[arXiv:hep-ph/0409040] . 

6. Q. H. Cao, R. Schwienhorst, J. A. Benitez, R. Brock and C. P. Yuan, Phys. Rev. D 72, 
094027 (2005) [arXiv:hep-ph/0504230] . 

7. V. M. Abazov et al . (DO Collaboration), Phys. Rev. Lett. 99, 191802 (2007) 
[arXiv:hep-ex/0702005] . 

8. V. M. Abazov et al. (DO Collaboration), accepted by Phys. Rev. Lett, larXiv: 0803^32561 
[hep-ex] . 

9. T. Tait and CP Yuan, Phys. Rev. D 63, 014018 (2001). 

10. V. M. Abazov et al. (DO Collaboration), accepted by Phys. Rev. D, arXiv: 0803 .07391 
[hep-ex] . 



